The Clusters AgeS Experiment (CASE). III. 
Analysis of the Eccentric Eclipsing Binary V32 in the Globular 

Cluster NGC 6397 1 

J. Kaluzny 2 , I. B. Thompson 3 , S. M. Rucinski 4 , W. Krzeminski 5 , 

ABSTRACT 

We present spectroscopic and photometric observations of the eclipsing binary 
V32 located in the central field of the globular cluster NGC 6397. The variable 
is a single-line spectroscopic binary with an orbital period of 9.8783 d and a 
large eccentricity of e = 0.32. Its systemic velocity (7 = 20.7 km s _1 ) and 
metallicity ([Fe/H] ~ -1.9) are both consistent with cluster membership. The 
primary component of the binary is located at the top of the main-sequence 
turn-off on the cluster color-magnitude diagram. Only a shallow primary eclipse 
is observed in the light curve. Based on stellar models for an age of 12 Gyr and the 
mass-function derived from the radial velocity curve, we estimate the masses to 
be M p = 0.79 M and M s = 0.23 M . The light curve of V32 can be reproduced 
by adopting R p = 1.569 R Q and R s = 0.236 R for the radii and i = 85.44 deg 
for the system inclination. The system geometry precludes observations of the 
secondary eclipse. The large eccentricity of the orbit is puzzling given that for 
metal poor, halo binaries the transition from circular to eccentric orbit occurs 
at an orbital period of about 20 days. We suppose that the orbit of V32 was 
modified relatively recently by dynamical interaction with other cluster star(s). 
An alternative explanation of the observed eccentricity calls for the presence of 
a third body in the system. 
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INTRODUCTION 



Little is known about the properties of main-sequence binary stars in globular clusters. 
These systems offer the potential for measurements of cluster age and distance independen t 
of cluster main-sequence fitting (see, for example, iPaczyhskil 119971 ; iThompson et al.ll200ll ). 
and as testbeds for theories of the stellar evolution of Population-II stars. 



The eclipsing binary V32-N6397 (hereafter V32) was discovered by lKaluzny et al.l (120061 ) 
during a survey for variable stars in the central field of the globular cluster NGC 6397. They 
observed only one shallow eclipse event with a depth of about 30 mmag. With V max = 16.13 
the variable is located at the very top of the cluster main sequence, and is a candidate 
member of NGC 6397. In this paper we report results of follow-up photometric and spectro- 
scopic observations of V32. Section [2] describes new photometry and a measurement of the 
spectroscopic orbit of V32. A self-consistent model of the system is presented in Section [3j 
Finally, Section H] briefly discusses the implications of the large measured eccentricity of the 
orbit of V32. 



2. SPECTROSCOPIC AND PHOTOMETRIC OBSERVATIONS 

Sp ectroscopic observatio ns of V32 stars were carried out with the MIKE echelle spectro- 



graph (IBernstein et al.ll2003l ) on the Magellan II (Clay) telescope at Las Campanas Obser- 
vatory. The data were collected during several observing runs between 2005 May and 2007 
August. For this analysis we use data obtained with the blue channel of MIKE covering 
the wavelength range 3350 A to 5000 A at a resolving power of A/AA ~ 38,000. All of 
the observations were obtained with a 0.7 x 5.0 arcsec slit and with 2x2 pixel binning. 
At 4380 A the resolution was ~ 2.7 pixels with a scale of 0.043 A/pixel. The spect ra were 



first p rocessed using a pipeline developed by Dan Kelson following the formalism of [Kelson 
( 2003 ) and then analyzed further using standard tasks in the IRAF/Echelle packag^. Each 



of the final individual spectra typically consisted of two 300 s exposures interlaced with an 
exposure of a thorium-argon lamp. We obtained 16 spectra of V32. 

Velocities were measured with the IRAF FXCOR package. For the temp late we used a 



single MIKE spectrum of the metal-poor subgiant HD 193901 ([Fe/H] = — 1.22. lTomkin et al. 



19921 ) with an adopted radial velocity of —172 km s The velocity measurements were 



made over the wavelength range 4000 A - 5000 A with the Balmer lines masked out of the 



6 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 
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template and object spectra. There was no evidence for a second velocity peak in any of 
the object spectra and we conclude that the variable is a single-line spectroscopic binary. 
The observations are presented in Table [T] which lists the heliocentric Julian Date (HJD) at 
mid-exposure, the velocities of the primary and errors of these velocities as returned by the 
FXCOR routine. 

The measured radial velocities of the primary were fitted with a non-linear least squares 
solution of the eccentric orbit using code written by G. Torres. The derived orbital elements 
are listed in the second column of Table [2J The quantity T\ is the moment of periastron 
passage and the remaining quantities have their standard meaning. The solution with 6 free 
parameters implies that a superior conjunction occured at HJD = 2453109.7977 ± 0.083. 
This is consistent with wit h the timing of the e clipse observed on the night of 2004 April 



11. Using photometry from iKafuzny et al.l (120061 ) supplemented with new observations (see 



below) we estimate that the observed eclipse was centered at HJD = 2453109.829. We then 
fitted the radial velocity curve by fixing the orbital period at the value which exactly predicts 
the moment of superior conjunction. By fitting 5 parameters we obtained P = 9.8779 d and 
T\ = 2453899.8248. This ephemeris predicts correctly the timing of the eclipse observed 
in the 2004 season. However, it leads to unacceptable phasing of some observations from 
the 2007 season which were obtained shortly before ingress to the primary eclipse. The 
shortest orbital period which gives a symmetric primary eclipse is P = 9.8783 d. For this 
period the observed moment of the observed superior conjunction can be reproduced for 
Ti = 2453899.8890. The solution with P and 7\ fixed at these values is listed in the third 
column of Table [2J The corresponding velocity curve and observed velocities are plotted in 

Fig.m 

At the beginning of the 2007 observing season we had in hand an approximate ephemeris 
for V32 and further photometric observations of a few predicted eclipses of the binary were 
obtained. The observations were made with the TEK5 camera on the 2.5m du Pont telescope 
at Las Campanas Obs ervatory. The i n strum ental setup and reduction methods were the same 



as these described in IKafuzny et al.l (120061 ). These observations were mostly hampered by 
poor weather but on one of the nights we obtained some useful data covering the phases 
immediately preceding the predicted ingress into primary eclipse. These observations helped 
to constrain the spectroscopic solution as described above. The l/-band light curve based on 
all available observations of V32 is shown in Fig. [2] for the region of the light curve near the 
observed eclipse. It is phased with the spectroscopic ephemeris listed in the third column 
of Table [2J The continuous line shows the synthetic light curve corresponding to the model 
presented in the next section. The full V-band observation data set is given in Table [3J 
The V/B — V color-magnitude diagram for NGC 6397 based on the du Pont photometric 
observations is presented in Fig. [3j 
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3. Properties of V32 



A direct determination of the absolute parameters of the components of V32 is hampered 
by the lack of velocity observations of the secondary. In the following analysis we assume 
that the variable is a member of the globular cluster NGC 6397. This conclusion is based on 
four arguments. First, the systemic velocit y of V32 (7 = 20.65 ±0.1 9 km s _1 ) agrees with the 
mean velocity of the cluster determined by lMeylan fe Mayorl (119911 ). V ra d = 18.1 ±0.1 km s _1 
with a central velocity dispersion of 4.5 km s _1 . Second, on the V/B — V diagram shown in 
Fig. [3] the binary is located at the top of the cluster turnoff, the primary is apparently just 
beginning to evolve onto the cluster sub-giant branch. Third, the variable is located on the 
sky only 51 arcsec from th e projected c enter of NGC 6397. The cluster itself has a half-light 
radius of = 140 arcsec ( Harris! 1996 ). Finally, we have estimated the metallicity of V32. 
We shifted the individual spectra of V32 to zero velocity, and averaged these shifted spectra 
to produce a mean spectrum for V32 . We measured the equivalent widths of 10 Fe I lines 
selected from iThompson et al.l (120081 ) in this mean spectrum, and compared these to the 



equiva lent widths measured in a grid of alpha-enhanced synthetic spectra from lCoelho et al. 
( 120051 ) linearly interpolated at T e ff = 6400 and log g = 4.3 for metallicities of [Fe/H] = -2.4, 
-2.0, -1.6, and -1.2. From this comparison we obtain [Fe/H] = -1.9 ± 0.1 (internal standard 
deviation) for V32. Recent measurements of the metallicity of NGC 6397 find [Fe/H] = -2.0 
( iGratton et al.ll2003l : iThevenin et al.ll200ll ). While our measurement is crude (a fine analysis 
of the mean spectrum is in preparation (Mc William et al. 2008)), the estimated metallicity 
of V32 is consistent with cluster membership. 

We adopt a cluster metallicity of [Fe/H] = —2.0 with an a-element enhancement of 
[ q/H] = +0.34. Rec ent estimates o f the age of NGC 639 7 span a range of 11.47 ± 0.47 Gyr 
(IHansen et al.l 120071 ) to 13.5 ±1.1 (IGratton et al.l (120031 ); a determination for models with 
di ffusion). We adopt a n age of 12.0 ±0.5 Gyr. Using the Dartmouth stellar model isochrones 
of botter et all J2007h for an age of 12.0 ± 0.5 Gyr, [Fe/H] = -2.0 and [a/H] = ±0.4, the 
location of the primary of V32 on the cluster V/ (B — V) diagram implies a primary mass of 
M p = 0.790 ± 0.009 M . Assuming for the moment an orbital inclination of i = 90 deg, we 
use the measured m ass function to derive a secondary mass of M s = 0.232 ± 0.003 M . The 
Potter et al.l (120071 ) models then suggest a secondary radius of R s = 0.236 R Q . 



We fit the light cu rve of V32 using the PHO EBE code dPrsa fe Zwitterl 120051) which is 
based on the models of Wilson fe Devinney (jl971 ) and Wilson ( 1979 ). Two free parameters 
were fit: the inclination i and the radius of the primary. The remaining parameters were fixed 
as follows: The orbital elements were adopted from the 3rd column of Table 2. The effective 
temperature of the primary was set at T p = 6490 K base d on an unr eddened (B — V ) p = 0.368 



and the empirical calibration of IWorthey fe Led (120061 ) . Following IGratton et al.l (120031 ) we 
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adopted E(B — V) = 0.183. The temperature of the secondary was set at T s = 3874 using 
the Dotter et al. (2007) models. The solution converged to i — 84.96 deg and R p = 1.53 R@. 
We then made one more iteration of the above procedure starting with i = 84.96 deg, M s = 
0.233 M Q and R s = 0.237 R Q . The solution converged to i — 85.03 deg and R p = 1.53 R Q . 
The calculated luminosity ratio in the l/-band is L p / L s = 664. In the middle of the primary 
eclipse the entire disk of the secondary is projected against disk of the primary - the eclipse 
is a transit. The secondary (occultation) eclipse is not observable. The estimated radius 
of the primary is consistent with its location above the main-sequence turnoff on the color- 
magnitude diagram of the cluster. For a mass of 0.790 M and an age of 12 Gyr the models 
of Dotter et al. (2007) predict a radius of r = 1.40 R Q . For the same age the radius reaches 
r = 1.54 R & for a mass of 0.793 M Q . Given the approximate nature of our analysis we 
consider the derived solution to be self-consistent. 

Clearly the limited data make a direct and accurate determination of the absolute 
parameters for components of V32 impossible. The light curve lacks a secondary eclipse 
and the photometric coverage of the primary eclipse is not complete. However, the deduced 
parameters of V32 reproduce well its observed radial velocity and light curves, and are 
self-consistent with evolutionary models of low mass stars. 



4. Discussion 



The properties of V32 are consistent with membership in the globular cluster NGC 6397. 
The system is composed of two main-sequence stars with an age of about 12 Gyr. The orbit 
shows a large eccentricity with e = 0.32. Such an eccentricity is unexpected is since tidal 
forces should circularize the orbit of the binary on a short time scale. The theoretical 
and observational aspects of the circularization of bin ary orbits for solar mass binaries are 
discussed in some det ail in iMeibom &: Mathieul (120051 ). They estimate, based on results of 
Latham et al. I ( 20021 ) . that halo binaries with periods shorter than 15.6^32 d should have 
circular orbits. For the old open clusters NGC 188 (age ~ 7 Gyr) and M67 (age ~ 4 Gyr) 
circularization has occured for periods shorter than 14.5 d and 12.1 d, respectively. 

Very little is known about globular cluster binaries with orbital periods longer than a 
few days. Several surveys for spectroscopic binaries among cluster giants were conducted in 
1980s and 1990s. They led to the detectio n of some candidates but generally spectroscopic 
orbits of these stars remain undetermined (jYan fc Cohen! Il996l . see also references therein). 
Four eclipsing binaries with periods 4 < P < 10 days wer e detected in the central part of 
47 Tuc in the HST/WFPC2 survey by lAlbrow et al.l (120011 ) using photometric observations. 
All of these have presumed circular orbits although observations do not cover the full orbital 
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cycle fo r two of the systems with the lo ngest periods. Extensive ground base sur veys of 
47 Tuc flWeldrake et alibopj ) and u Cen fcaluznv et aDhood ; IWeldrake et"alll2007h led to 
the detection of two eclipsing binaries with orbital periods exceeding 10 days. Follow up 
observations conducted by our group show that 47 Tuc-V69 and u Cen-V406 have eccentric 
orbits with P = 29 d and P = 71 d, respectively (Thompson et al. 2008, in preparation). 
All of the other cluster eclipsing binaries reported so far have orbital periods less than 10 
days and show circular orbits. 

The eccentric orbit of V32 can possibly be due to a relatively recent dy namical interac- 
tion o f the binary with other cluster star(s). NGC 6397 has a "collapsed" core (jDjorgovski fc King 
19861 ) containing a population of about 20 X-ray sources including 9 candidat e cataclysmic 



varia bl es, a millisecon d pulsa r and several candidate BY Dra-type close binaries (IGrindlay et al. 



200ll ). iKaluzny et al.l (120061 ) reported the detection of 9 eclipsing binaries and 6 candidates 
for ellipsoidal binaries located in the central part of the cluster. A dynamical interaction 
capable of transforming the circular orbit of V32 into a highly eccentric one would also likely 
significantly disturb the systemic velocity of the binary. Our data show that the radial ve- 
locity of V32 is not unusual for an object from the central part of the cluster. Unfortunately, 
the transverse velocity is V32 is unknown as the binary i s not included in the recent proper 
motion survey of NGC 6397 conducted by iKalirai et al.l (120071 ). 



As an alternative explanation of the eccentric orbit of V32 we consider the possibility 
that the system is in fact a hierarchical triple. The third component located on the outer 
orbit is capable of generating eccentricity in the inner binar y. Exte nsive reviews of the 



evolut ion of the orbits of binary stars have been presented by lEggletonl (120061 ) and iMazeh 
( 120081 ). Numerical simulations of stellar clusters pr edict dynamical formation of triple stars, 
espec ially in the presence of primordial binaries ([McMillan et al.l Il990l ; iHeggie fc Aarseth 
19921 ). Ex amples of known trip le stars in globular clusters include the pulsar PSR B1620- 
26 in M4 (IThorsett et al.lll993l ) and possibly the ultracompact X-ray binary 4U 1820-303 
dzdziarski et al. Il2007h . The hypothesis that V32 is a triple can be tested by spectroscopic 
observations aimed at detection of variability of the systemic velocity of the putative inner 
binary. We see no systemic velocity residuals from the orbit give in Table 3 over the 2.2 year 
duration of our velocity observations. 

Finally we comment on possibility that the ec centric orbit of V32 i s rela ted to its rather 
low mass ratio q = M s /M p rj 0.30. According to iMathieu fc Mazehl (119881 ) the time scale 



of circularization of a binary is related to q by the relation: r circ oc q 2 ' 3 (l + g -1 ) 5 ' 3 . This 
corresponds to an increase of r C j rc of a factor of 1.6 as the mass ratio ranges from q = 1 to 
q = 0.3. For low val ues of q the contribution of the secondary to the circularization process 
becomes negligible (lMazehll2008l ). At the same time T C j rc is a strong function of the orbital 
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period P, with r C j rc oc P 16 / 3 ( Zahnlll9771 ). The net result is that the circularization period 
will lengthen as the mass ratio decreases from for q = 1 and q = 0.3, especially for long 
period syste ms. It is appropriate to note at this point that that the sample of halo stars 
presented by lLatham et al. I (120021 ) consists entirely of single line binaries and hence objects 
with mass ratios noticeably lower than unity. More globular cluster binaries with period 
of the order of 10 days have to be detected and analyzed before V32 can be called a truly 
unusual system. The CASE group is in a process of collecting data which can shed more 
light on the relation between eccentricity and orbital period for binaries in globular clusters. 



JK and WK were supported by the grant 1 P03D 001 28 from the Ministry of Science and 
Higher Education, Poland. Research of JK is also supported by the Foundation for Polish 
Science through the grant MISTRZ. IBT was supported by NSF grant AST-0507325. Sup- 
port from the Natural Sciences and Engineering Council of Canada to SMR is acknowledged 
with gratitude. It is a pleasure to thank Willy Torres for sharing his code with us. 



REFERENCES 

Albrow, M.D., Gilliland, R.L., Brown, T.M., Edmonds, P.D., Guhathakurta, P., & Saraje- 
dini, A 2001, ApJ, 559, 1060 

Bernstein, R., Shectman, S. A., Gunnels, S. M., Mochnacki, S., & Athey, A. E. 2003, Instru- 
ment Design and Performance for Optical/Infrared Ground-based Telescopes. Edited 
by lye, Masanori; Moorwood, Alan F. M. Proceedings of the SPIE, 4841, 1694 

Coelho, P., Barbuy, B., Melendez, J, Schiavon, R.P, & Castilho, B.V. 2005, A&A, 443, 735 

Djorgovski, S., & King, I.R. 1986, ApJ, 305, L61 

Dotter, A., Chaboyer, B., Jevremovi, D., Baron, E., Ferguson, J. W., Sarajedini, A., Ander- 
son, J. 2007, AJ, 134, 376 

Eggleton, P. 2006, Evolutionary Processes in Binary and Multiple Stars, Cambridge Astro- 
physics Series, No. 40 

Gratton, R.G., Bragaglia, A., Carretta, E., Clementini, G., Desidera, S., Grundahl, F., & 
Lucatello, S. 2003, A&A, 408, 529 

Grindlay, J.E., Heinke, CO., Edmonds, P.D., Murray, S.S., & Cool, A.M. 2001, ApJ, 563, 
53 



- 8- 



Harris, W. E. 1996, AJ, 112, 1487 

Heggie, D.C., & Aarseth, S.J. 1992, MNRAS, 257, 513 

Hansen, B. M. S., et al. 2007, ApJ, 671, 380 

Kaluzny, J., Olech, A., Thompson, LB., Pych, W., Krzeminski, W., Schwarzenberg-Czerny, 
A. 2004, A&A, 424, 1101 

Kaluzny, J., Thompson, LB., Krzeminski, W., Schwarzenberg-Czerny, A. 2006, MNRAS, 
365, 548 

Kalirai, J.S., et al. 2007, ApJ, 657, L93 
Kelson D. D. 2003, PASP, 115, 688 
Latham, D.W. et al. 2002,AJ, 124, 1144 
Mathieu, R.D., & Mazeh, T. 1988, ApJ, 326, 256 

Mazeh, T. 2008, in Tidal Effect in Stars, Planets, and Disks, Ed. M.-J. Goupil and J.-P. 
Zahn, EAS Publications Series, in press, preprint larXiv:0801 .0134 

McMillan, S, Hutt, P, & Makino, J. 1990, ApJ, 362, 522 

Meibom, S., Mathieu, R.D. 2005, ApJ, 620, 970 

Meylan, G., & Mayor, M. 1991, A&A,250, 113 

Paczyhski, B. 1997, in The Extragalactic Distance Scale STScI Symposium (Cambridge: 
Cambridge University Press) 

Prsa, A., Zwitter, T. 2005, ApJ, 628, 426 

Thevenin, F., Charbonnel, C, de Freitas Pacheco, J. A., Idiart, T. P., Jasniewicz, G., de 
Laverny, P., & Plez, B. 2001, A&A, 373, 905 

Thompson, I. B., Kaluzny, J., Pych, W., Burley, G., Krzeminski, W., Paczyhski, B., Persson, 
S. E., & Preston, G. W., AJ, 121, 3089 

Thompson, I. B., Ivans, I. I., Bisterzo, S., Sneden, C., Gallino, R., Vauclair, S., Burley, G. 
S., Shectma, S. A., & Preston, G. W. 2008, ApJ, in press, preprint arXiv:0712:3228 

Tomkin, J., Lemke, M., Lambert, D.L., & Sneden, C. 1992, AJ, 104, 1568 



- 9- 

Thorsett, S.E., Arzoumanian, Z., & Taylor, J.H. 1993, ApJ, 412, L33 

Weldrake, D.T.F., Sackett, P.D., Bridges, T.J., & Freeman, K.C. 2004, AJ, 128, 736 

Weldrake, D.T.F., Sackett, P.D., & Bridges, T.J. 2007, AJ, 133, 1447 

Wilson, R. E., & Devinney, E. J. 1971, ApJ, 166, 605 

Wilson, R. E. 1979, ApJ, 234, 1054 

Worthey, G., & Lee, H.-C. 2006, preprint dastro-ph/0604590| 
Yan, L., & Cohen, J.G. 1996, AJ, 112, 1489 
Zahn, J.-R, 1977, A&A, 57, 383 

Zdziarski, A.A., Wen, L., & Gierliofiski, M. 2007, MNRAS, 377, 1006 
Zucker, S. & Mazeh, T. 1994, ApJ, 420, 806 



This preprint was prepared with the AAS IATgX macros v5.0. 



Table 1. Velocity Observations of V32 
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Table 2. Orbital Parameters for V32 



Parameter Value Value 

P (days) 9.8783 ± 0.0011 9.8783 a 

Ti (HJD-245 0000) 3899.832 ± 0.083 3899.889 a 

7 (km s- 1 ) 20.67 ± 0.22 20.71 ± 0.18 

K p (km s' 1 ) 23.92 ± 0.30 24.00 ± 0.25 

e 0.320 ± 0.011 0.322 ± 0.009 

uj (deg) 72.7 ±3.0 74.74 ± 0.7 

(O - C) rms (km s' 1 ) 0.60 0.62 
Derived quantities: 

f(M) sin 3 i (M & ) 0.01192 ± 0.00046 0.01201 ± 0.00042 

A p sin(i) (R Q ) 3.079 ± 0.039 3.086 ± 0.036 



a Fixed based on photometry. 
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Table 3. V-band Photometric Observations of V32 
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Note. — This Table is pub- 
lished in its entirety in the elec- 
tronic addition of the Astro- 
nomical Journal. A portion 
here is shown for guidance re- 
garding its form and content. 
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Fig. 1.— The velocity curve for orbital solution for V32, together with the individual 
observed velocities. The O-C residuals from the solution are shown in the lower panel. 
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Fig. 2. — The phased observed V light curve of V32. The synthetic curve corresponding to 
the model described in Sec. 3 is shown in the upper panel. 
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Fig. 3.— Position of V32 in the V/B - V diagram for NGC 6397. 



